The effect of temperature jumps during polymer crystallization 
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Temperature changes during the growth of lamellar polymer crystals give rise to steps on the 
surface of the crystals. It has recently been suggested that these steps could provide important 
insights into the mechanism of polymer crystallization. In particular, a characterization of the 
profiles of these steps might reveal the fixed-point attractor that underlies a recently proposed 
crystallization mechanism. Here we examine this hypothesis by performing simulations of such 
temperature jumps using the Sadler-Gilmer model. We find that for this model the step profiles do 
reveal the fixed-point attractor. However, for temperature decreases they also reflect the rounding 
of the crystal edge that occurs in this model and for temperature increases they also reflect the 
fluctuations in the thickness present in the crystal. We discuss the implications of these results for 
the interpretation of experimental step profiles. 
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I. INTRODUCTION 

Upon crystallization from solution and the melt many 
polymers form lamellae where the polymer chain tra- 
verses the thin dimension of the crystal many times fold- 
ing back on itself at each surfacelm (The crystal ge- 
ometry is shown by the example configuration in Fig- 
ure [l]) Although lamellar crystals were first observed 
over forty years ago their physical origin is still con- 
troversial. In particular the explanations for the depen- 
dence of the lamellar thickness on temperature offered by 
the two dominant theoretical approaches — the Lauritzen- 



and the entropic 
— differ greatly.EJ 



Hoffman surface nucleation theory . 
barrier model of Sadler and Gilmeit 

One of the common features of the two theories is 
that they both argue that the observed crystal thick- 
ness is close to the thickness at which the crystal growth 
rate is a maximum. However, recently a new descrip- 
tion of thc^ mechanism of thickness selection has been 
presented. IlirEj In this approach the observed thickness 
corresponds instead to the one thickness, I**, at which 
growth with constant thickness can occur. Crystals ini- 
tially thicker (thinner) than I** will thin (thicken) as the 
crystals grow until the thickness I** is reached. This dy- 
namical convergence can be described by a fixed-point 
attractor which relates the thickness of a layer to the 
thickness of the previous layer. The value of the thick- 
ness at the fixed point is £**. 

This mechanism has been feujad for two simple mod- 
els of polymer crystallizationJilriJ In the first model the 
polymer crystal grows, as in the Lauritzen- Hoffman the- 
ory, by the successive deposition of stems (a stem is a 
straight portion of the polymer chain that traverses-tha 
thin dimension of the lamella) across the growth faceJiMa 
A configuration produced by this model is shown in Fig- 
ure [I] to illustrate the mechanism. The crystal thins down 
from the initial thickness to I** within five to ten layers 



and then continues to grow at that thickness The second 
model is that used by Sadler and Gilmer, the behaviour 
of which they interpreted in terms of an entropic barrier. 
In this Sadler-Gilmer (SG) model the connectivity of the 
polymer is modelled implicitly, the growth face can be 
rough, and lateral correlations along the growth face can 
be weak. That we find the same mechanism in these two 
very different models is a sign of its generality. 



direction of 
crystal growth 




FIG. 1. Cut through a polymer crystal which was pro- 
duced by the growth of twenty successive layers on a sur- 
face with a uniform thickness of 50 units using the model 
described in Refs. and Folds occur at the top and 
bottom surfaces, and the stems are represented by vertical 
cuboids. The cut is 16 stems wide. 



Further support for the new mechanism is provided by 
the experimental observation that a temperature chpasft 
during crystallization produces a step on a lamellaliirEj 
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This step is a result of the thickness of the crystal dynam- 
ically converging to I** for the new temperature as the 
crystal grows. Furthermore it has been suggested that a 
detailed characterization of the step profiles by atomic- 
force microscopy could allow the fixed-point, attractors 
that underlie the mechanism to be obtainedJijtLj In this 
paper we examine this suggestion more carefully by per- 
forming simulations of temperature jumps for the Sadler- 
Gilmer model.li3 In particular, we investigate the effect 
that rounding of the crystal profile near to the growth 
face and fluctuations in the crystal thickness may have 
on the shape of the steps. It is hoped that this work will 
aid the experimental interpretation of step profiles. 




FIG. 2. A schematic picture of a two-dimensional slice 
(perpendicular to the growth face) through a lamellar poly- 
mer crystal which forms the basis of the two-dimensional ver- 
sion of the Sadler-Gilmer model. The three possible changes 
in configuration allowed by the model are shown (the dashed 
lines represent the outline of the possible new configura- 
tions). 

II. METHODS 

In the SG model the growth of a polymer crystal re- 
sults from the attachment and detachment of polymer 
units at the growth face. The rules that govern the sites 
at which these processes can occur are designed to mimic 
the effects of the chain connectivity. In the original three- 
dimensional version of the model, kinetic Monte Carlo 
simulations were performed to obtain many realizations 
of the polymer crystals that result. Averages were then 
taken over_these configurations to get the properties of 
the modeLQ Under many conditions the growth face is 
rough and the correlations between stems in the direction 
parallel to the growth face are weak. Therefore, an even 



simpler two-dimensional version of the model was devel- 
oped in which lateral correlations are neglected entirely, 
and only a slice through the polymer flpystal perpendic- 
ular to the growth face is considered flu The behaviour 
of this new model was found to be very similar to the 
original three-dimensional model. 

The geometry of the model is shown in Figure ||. 
Changes in configuration can only occur at the outer- 
most stem and stems behind the growth face are 'pinned' 
because of the chain connectivity. There are three ways 
that a polymer unit can be added to or removed from the 
crystal: (1) The outermost stem can increase in length 
upwards. (2) A new stem can be initiated at the base of 
the previous stem. (3) A polymer unit can be removed 
from the top of the outermost stem. 

The ratio of the rate constants for attachment (fc + ) 
and detachment (k~) of a polymer unit are related to 
the thermodynamics of the model through 

k + /k~ = exp(-AF/fcT), (I) 

where AF is the change in free energy on addition of a 
particular polymer unit. The above equation only defines 
the relative rates and not how the the free energy change 
is apportioned between the forward and backward rate 
constants. We follow Sadler and Gilmer and choose k + 
to be constant. We use l/k + as our unit of time. 

In the model the energy of interaction between two 
adjacent crystal units is -e and the change in entropy 
on melting of the crystal is given by AS = AH/T m = 
2e/T m , where T m is the melting temperature (of an in- 
finitely thick crystal) and AH is the change in enthalpy. 
It is assumed that AS is independent of temperature. 
Here, as with Sadler and Gilmer, we do not include any 
contribution from chain folds to the thermodynamics. 

From the above considerations it follows that the rate 
constants for detachment of polymer units are given by 

k~(i, j) = k + exp(2e/kT m - 2e/kT) i^l,i<j (2) 
k~(i, j) = k + cxp(2e/kT m - e/kT) i = l,i>j, (3) 

where i is the length of the outermost stem and j the 
length of the stem in the previous layer. The first term 
in the exponents is due to the gain in entropy as a result 
of the removal of a unit from the crystal, and the second 
term is due to the loss of contacts between the removed 
unit and the rest of the crystal. 

There are two ways to examine the behaviour of the 
model. In one approach the model is formulated in terms 
of a set of rate equations which can easily be solved 
numerically to yield the steady-state solution of the 
modelfltj This is the method that we used for the most 
part in our previous study of the Sadler-Gilmer modelxiH 
However, as we wish to examine the evolution of the sys- 
tem towards the steady state, we use kinetic Monte Carlo 
to grow a set of representative crystals.EJ In this work 
we deliberately start growing these crystals from a well- 
defined non-steady-state initial configuration — either a 
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crystal of constant thickness different from I**, or a crys- 
tal grown at a different temperature. Averages are then 
taken over these crystals to obtain information about the 
convergence of the system towards the steady state. 

At each step in the kinetic Monte Carlo simulation a 
state, b, is randomly chosen from the three states con- 
nected to the current state, a, with a probability given 
by 

P ab = (4) 

l^b' K ab' 

and update the time by an increment 



where p is a random number in the range [0,1]. Depend- 
ing on the conditions we use from tens of thousands to 
millions of steps to grow each individual crystal, and then 
take averages over many thousands of crystals. 

The version of the SG model used here has two vari- 
ables: kT m /e and-Ij/Tm. Here, as in our previous work 
on the SG modeljlj we use kT m /e = 0.5, unless other- 
wise stated. Sadler and Gilmer have shown that the basic 
properties of the model were independent of the value-pf 
kT m /e within the parameter range that they studied. □ 

III. RESULTS 

Before considering simulations of actual temperature 
jumps we examine a slightly simpler case, namely the 
convergence of the crystal thickness to I** when the initial 
configuration is a crystal of constant thickness different 
from I**. These cases will provide a useful comparison 
to the convergence to a new I** caused by a change in 
temperature. 

Figure [}] shows two example crystal profiles when the 
the thickness of the initial crystal is larger than I**. It 
can be clearly seen that the crystal quickly converges to 
I** whatever the initial thickness of the crystal leading 
to a downward step on the surface. For the crystal that 
is initially twenty units thick there is little backward mo- 
tion of the growth front in the early stages of the simu- 
lations because a thick crystal such as this is very stable. 
Therefore, after growth only the outer layer of the ini- 
tial crystal has a thickness different from its initial value. 
The step has a sharp downward edge. 

By contrast, the step on the crystal that was initially 
ten units thick is more rounded (Figure |b). The outer 
layers of the initial crystal are now noticeably less than 
ten up to about five layers back into the initial crystal 
and the curvature of the step changes from negative to 
positive around the initial position of the edge of the 
crystal. 

It is also worth noting that we previously found that 
the initial growth rate increases with the thickness of the 



initial crystal.113 The cause of this behaviour is similar 
to that for the more pronounced rounding of the steps 
that result from growth on thinner crystals. Both are 
related to the smaller number of detachment steps in the 
kinetic Monte Carlo simulations when the initial crystal 
is thicker because of the greater thermodynamic driving 
force for growth. 
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FIG. 3. Average crystal profiles that result from the 
growth of an initial crystal that has a constant thickness of 
(a) 20, (b) 10 and (c) 7 polymer units. For (a) and the solid 
line in (b) the zero of position is defined as the edge of the 
initial crystal. For (c) and the dashed line in (b) the zero of 
position is defined as the minimum position of the growth 
front during the simulation. T — 0.95T m . 
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FIG. 4. Dependence of the thickness of a layer on the 
thickness of the previous layer. The lines with data points 
in (a)-fc), lj), are obtained from the crystal profiles in 

Figure |3|(a)-(c). The arrows on these lines mark the direc- 
tion of growth and the label '0' marks the point at which 
j=0. The dotted lines are y = x and l n (l n -i), as labelled. 
T = 0.95T m . 

In Figure |j|b, as well as the average profile with respect 
to the position of the edge of the initial crystal, we also 
show the average profile measured with respect to the 
minimum position of the growth front in each individ- 
ual crystal. Now the sharp edge of the step is regained. 



Therefore, the rounding of the step for the space-fixed 
profile can be understood to result from the variation in 
the amount by which the growth face transiently retreats 
back into the initial crystal. 

To use the step profile to reveal the fixed-point at- 
tractor that describes the convergence to I** we plot in 
Figure |^ the thickness of a layer against the thickness of 
the previous layer (i.e. the points lj) where j is the 

position of a layer) as we go through the step. The initial 
points from layers behind the step are all at (k n it, h n it)> 
where l- m i t is the thickness of the initial crystal. As one 
passes through the step the points leave (Zinit, 'init) and 
follow a path which converges to (/** , I**). It is the nature 
of this path that is of interest. In particular, the assump- 
tion of our previous suggestion that the steps produced 
by temperature changes could provide insight into the 
mechanisms of polymer crystallization was that this path 
should follow the fixed-point attractor, i n (Z„_x). l n (l n — \) 
is defined as the average thickness of a layer n in the 
bulk of the crystal given that the thickness of the pre- 
vious layer is l n -i and is obtained from the steady-state 
solution of the rate equations describing the SG model 
(for details see Ref. [b|). 

In Figure ^ we also plot l n (l n -i) to compare with the 
path taken by (lj-x,lj). For the crystal that is initially 
twenty units thick the path jumps from (20,20) to the 
fixed-point attractor in two steps and then follows it 
down to the fixed point. The intermediate point is due to 
the slight rounding of the outer layer of the initial crystal. 
The plot for the crystal that is initially 10 units thick is 
similar, except that the number of steps taken to reach 
the fixed-point attractor is larger because of the greater 
rounding of the the initial crystal. However, when the 
position is measured with respect to the minimum po- 
sition of the growth front the path of immedi- 
ately steps onto the fixed-point attractor from the point 
(10,10). 

The third step profile in Figure |^ shows the profile 
when the thickness of the initial crystal is less than I**. 
In this case the initial crystal is unstable with respect to 
the melt/solution and so the growth face retreats (Figure 
||a). Only once a fluctuation to a greater thickness oc- 
curs does the growth face begin to advance. Generation 
of this fluctuation is a slow process because there is an 
energetic cost associated with the stems which overhang 
the previous layer. As the distance that the growth face 
initially goes backwards has a wide variation between in- 
dividual simulations, the gradient of the step in a space- 
fixed frame of reference is very shallow. Only when the 
profile is measured with respect to the minimum position 
of the growth face does a clear picture of the thickening 
emerge. In this frame of reference there is a sharp up- 
ward step and the thickness quickly reaches I** (Figure 
^c). Furthermore, the path of (lj-i,lj) jumps straight 
from (7,7) onto the thickening branch of the fixed-point 
attractor (Figure ||c) . 

At this point it is right to consider which frame of 
reference — space-fixed or fixed with respect to the min- 
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imum position of the growth face — is more appropriate 
to the step profiles on real polymer crystals. The profiles 
resulting from the two frames of reference represent two 
limits in the degree of correlations between events along 
the growth face. The space- fixed frame of reference main- 
tains the assumption of the two-dimensional SG model 
that there are no correlations between adjacent stems 
along the growth face. It is this assumption that leads 
to the variation in the distance by which the growth face 
retreats. However, this assumption is not always a good 
one. In particular, it seems likely that the nucleation of 
a thicker region would propagate laterally. Therefore, we 
expect the degree of correlation to be closer to the limit 
obtained from using a frame of reference fixed with re- 
spect to the minimum position of the growth face. This 
latter approach is equivalent to assuming that the line 
of the step along the fold surface is straight, rather than 
rough. 



(a) 




50000 100000 150000 200000 250000 300000 

k*t 

FIG. 5. Position of the crystal growth front as a function 
of time. In (a) the two lines are for growth from an initial 
crystal of thickness 7 at T = 0.95 T m and for growth after 
a temperature jump from 0.935 T m to 0.95 T m , as labelled, 
(b) shows the growth after a temperature jump from 0.9 T m 
to 0.95 T m . 

Having shown that for the situation when growth oc- 



curs from a crystal that is of constant thickness a plot of 
lj) can reveal the fixed-point attractor that under- 
lies the growth mechanism of polymer crystals, we now 
proceed to consider the steps that result from temper- 
ature jumps (from a temperature T± to a temperature 
T2) during growth. As the thickness of polymer crystals 
are approximately inversely proportional to the degree of 
supercoolingJ13 a decrease in temperature will lead to a 
downward step, and an increase in temperature will lead 
to an upward step. 
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FIG. 6. Average crystal profiles that result from tem- 
perature jumps (a) from 0.95 T m to 0.9 T m with fcT m /e=0.5 
and 0.2, as labelled, and (b) from 0.9 T m (solid line with di- 
amonds), 0.935 T m (dashed line with crosses) and 0.945 T m 
(dotted line with squares) to 0.95 T m . The zero of position 
is defined in (a) as the position of the growth face when the 
temperature jump occurs and in (b) as the minimum posi- 
tion of the growth face during the simulation. In (a) both 
the crystal profile before and after the temperature jump are 
shown. 
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FIG. 7. Dependence of the thickness of a layer on the 
thickness of the previous layer. The lines with data points, 
lj) are obtained from the step profiles in Figure pi The 
plots correspond to a temperature jump from 0.95 T m to 
0.9 T m with (a) fcT m /e=0.5 and (b) fcT m /e=0.2 and (c) to 
temperature jumps from 0.9 T m (solid line with diamonds), 
0.935 Tm (dashed line with crosses) and 0.945 T m (dotted line 
with squares) to 0.95 T m . The arrows on these lines mark the 
direction of growth and the label '0' marks the point at which 
j=0. In (a) and (b) the dashed line is for the profile of the 
crystal edge before the temperature jump and the solid line 
for the step profile after the temperature jump. The dotted 
lines are y = x, ln(ln-i,T2) and l' n (l n -\,T\), as labelled. 



We consider the effects of a decrease in temperature 
first, fn Figure ||a we show the average crystal profile 
before and after the decrease in temperature. Firstly, 
unlike the situation considered above, the new growth 
after the temperature change is on a crystal where there 
are variations in the layer thickness and where the lay- 
ers near to the growth face are on average thinner. The 
rounded profile at the growing edge of the crystal is a 
characteristic property of the SG model and plays a key 
role in Sadler and Gilmer's explanation ofjDjOlymer crys- 
tallization in terms of an cntropic barrier .013 

In the growth after the temperature jump not all of 
the rounding present at the edge of the crystal at the 
time of the jump is removed. The profile of the result- 
ing step initially follows the profile of the rounded edge, 
before changing curvature and smoothly converging to 
I** for the new temperature (Figure ^a). The path of 
for the step initially leaves (P*(Ti), Z**(Ti)) and 
follows the same path as (lj-i, lj) for the crystal edge at 
T\. Only when this line meets Z„(Z„_i,T 2 ) does the path 
change slope and follow the fixed-point attractor down 
to (l**(T 2 ),l**(T 2 )) (Figure 0a). 

This basic scenario holds for all temperature decreases. 
The main differences are only in the degree to which the 
step reflects the rounding of the crystal edge, which in 
turn depends upon the relative slopes of l n (l n -i, T 2 ) and 
the path of (lj-i,lj) for the crystal edge at T\. For 
instance, if the product of the slopes is one then the 
crossover will occur halfway between l**(Ti) and l**(T 2 ), 
and for the example shown in Figure ^a this is approxi- 
mately the case. 

The parameters in the model that can affect the two 
slopes are kT m /e, T\ and T 2 . We do not intend to survey 
the full parameter space, but instead just comment on 
the effect of varying each parameter alone on the exam- 
ple in Figure |7]a. For example, if we decrease kT m /e the 
slope of the fixed-point attractor becomes closer to one. 
Therefore, the convergence of the thickness to l**(T 2 ) is 
more gradual (Figure ||a) and the path of (lj-±,lj) fol- 
lows-the fixed point attractor to a greater extent (Figure 
^b) .Ej However, increasing T\ has an opposite effect. It 
makes the slope of (lj-i, lj) for the crystal edge closer to 
one — the rounding is more gradual and extends deeper 
into the crystal away from the crystal edge. Therefore, 
the path of (Zj_i,Zj) follows the fixed point attractor to 
a lesser extent. Finally, changing T 2 has relatively little 
effect on the relative slopes and so the crossover remains 
roughly midway between 1**{T\) and l**(T 2 ). 

In Figure ^b we show three examples of steps that re- 
sult from temperature increases. In these cases a crystal 
of thickness of l**(T±) is unstable with respect to the 
melt /solution at T 2 and so the crystal growth face ini- 
tially retreats after the temperature jump (Figure ||). 
In one of the cases (Ti = 0.935) we chose T\ so that 
l**(Ti) ss 7 enabling us to make a comparison with the 
step that was produced when the initial crystal had a 
constant thickness of 7 units. From Figure ||a we can see 
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that growth begins markedly earlier when there is a tem- 
perature jump. The reason for this difference becomes 
clear when we examine the step profiles shown in Figure 
||b. The majority of the step is behind the minimum posi- 
tion of the growth face after the temperature jump. The 
growth face must retreat until it reaches a region of the 
crystal where the layer thickness is larger than the aver- 
age value at T\. New growth then begins from this posi- 
tion and the crystal thickens the small amount necessary 
to reach l**(T 2 ). Only for this new growth does 
follow the fixed-point attractor (Figure |c). Growth is 
more rapid than for the case where growth is from an 
initial crystal of constant thickness because encounter- 
ing already present fluctuations during the retreat of the 
growth face is a more common event than the generation 
of new fluctuations to larger thickness. 

Interestingly, the part of the step resulting from al- 
ready present fluctuations has a well-defined behaviour. 
In the bulk of the crystal (where the influence of the 
growth face can no longer be felt) there is a symmetry 
between the directions towards and away from the growth 
face. Therefore, i n _ i(l n ) — l n (l n -i)', i.e. in the bulk of 
the crystal the dependence of the thickness of layer n—1 
on the thickness of layer n is the same as the dependence 
of the thickness of layer n on the thickness of layer n—1. 
Given that at the minimum position of the growth face 
there is a fluctuation to a certain amount larger than 
I**, the layers behind this would therefore be expected 
to obey i n _i(/ n ,Ti). So, it is unsurprising that a plot 
of for j < follows the fixed point attractor for 

T\. From this it simply follows that a plot of lj) for 
j < follows the curve, l' n (l n -i,Ti), formed by reflect- 
ing Z„(Z n _i,Ti) in y = x (Figure 0c). For j > (lj-i,lj) 
jumps from i^(/ n _i, T\) onto the fixed point attractor for 
T 2 . 

From Figures ||, ||b and 0c we can ascertain some of the 
effects of changing the size of the temperature increase. 
When the temperature jump is larger the growth face 
retreats further because a larger fluctuation in the thick- 
ness needs to be encountered to nucleate new growth. 
Furthermore, the larger temperature jump the more of 
the step is a result of new growth. Interestingly, for the 
smallest temperature increase the growth face retreats to 
a fluctuation in thickness which is on average larger than 
Z**(T 2 ). Therefore, the profile displays a lip (Figure ||b) 
and in the new growth the crystal thins slightly to reach 
the new I**. 

The effect of kT m /e is somewhat similar to the effect 
of the magnitude of the temperature increase. A larger 
proportion of the step is the result of new growth when 
kT m /e is smaller. Presumably, this is because there is 
less variation in the stem length for smaller values of 
kT m /eB 



IV. DISCUSSION 

In this paper we have seen that the step profiles ob- 
tained by simulations of temperature jumps using the SG 
model do reveal information about the fixed-point attrac- 
tor that we have recently argued underlies the mechanism 
of polymer crystallization.EJTl3 This strengthens our sug- 
gestion that temperature jump experiments can provide 
insight into the physics of polymer crystallization. How- 
ever, the step profiles for temperature decreases also re- 
flect the rounding of the crystal edge and those for tem- 
perature increases also reflect the variations in the stem 
length that are present in the crystals. These additional 
features mean that the interpretation of experimental 
step profiles in terms of fixed-point attractor curves re- 
quires considerable care. This disadvantage is partly off- 
set by the the fact that an experimental study of the 
step profiles could reveal information about aspects of 
polymer crystallization not originally anticipated. 

However, we note that it is not clear to what extent 
these additional features occur for the steps on real poly- 
mer crystals — they may just reflect some of the simpli- 
fications in the SG model. Firstly, although the crystal 
edge is always rounded in the SG model, there is, as far as 
we are aware, not yet any direct experimental evidence 
of this effect occurring in real polymer crystals under 
normal crystallization conditions. Secondly, in an alter- 
native model of polymer crys talli zation, rounding only 
occurs at small supercoolingsJliO 

Secondly, it may be that the SG model overestimates 
the roughness of the fold surface of polymer crystals. In 
the SG model changes in the length of a stem can oc- 
cur only when that stem is at the growth face. Once 
the growth face has passed through a region the stem 
length variations that arise from the kinetics of growth 
are frozen in. However, it may be that in real crystals 
there are annealing mechanisms (which only require lo- 
cal motion of stems) that can act to reduce the mag- 
nitude of these fluctuations to their equilibrium values. 
Such processes might make the crystal growth following 
a temperature increase more similar to the growth from 
an initial crystal of constant thickness less than Z**.EJ 
Recent atomic-force microscopy experiments on the fold 
surface of polyethyleneEJa shed some light on this issue. 
Only for a minority of crystal could the folds be clearly 
resolved. It was suggested that the lack of clear order for 
the majority of crystals resulted from differences in the 
heights of the folds, i.e. there is some surface roughness. 
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